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Abstract 

We calculate the low-energy parameter S in QCD, which is also known 
as Lio, and the pion decay constant f^^ using inhomogeneous Bethe-Salpeter 
equation in improved ladder approximation. To extract these quantities we 
calculate the — ^" two-point function, Iivv{(f') — ^aa{<i'^), in space-like 
region. We obtain S = 0.43 ~ 0.48, which is about 30% larger than the 
experimental value. The calculated fj^ is well consistent with the result by 
solving the homogeneous Bethe-Salpeter equation for pion. We also evaluate 
S parameter in SU (3) gauge theory with doublets of fermions in connec- 
tion with walking technicolor model, and find that the value of S/No hardly 
depends on Nr,. 
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1 Introduction 



There is much interest to investigate the low-energy dynamics of QCD. We can see 
its property from the low-energy parameters of the effective Lagrangian such as Li, 
L2, Lio, introduced by Gasser-Leutwyler[l|], as well as the pion decay constant 
The parameter Lio is related to the S parameter, which expresses one of the 
oblique corrections p|, ^, ^ in electroweak theory. 

These low-energy parameters are calculated by various models. For example the 
free quark model gives a half of the experimental value for QCD S parameter. It is 
shown that the parameter Lio as well as the other parameters (Li, L2 and Lg) are 
saturated by the contribution from the low-lying vector and axial-vector mesons, p 
and ai.0 In Ref. 0, |^, based on the nonlocal constituent-quark model, QCD S 
parameter is calculated using a momentum dependent quark mass function. The 
estimated QCD 5* parameter well reproduces the experimental value. In Ref. 
they argue the corrections to the free quark loop diagram, and conclude that it is 
important to include the interaction which forms bound states as well as corrections 
to the quark propagator. 

In this paper, we calculate QCD 5* parameter (i.e., Lio) using the inhomogeneous 
Bethe-Salpeter (BS) equation in the improved ladder approximation. The quark 
mass function is consistently calculated by Schwinger- Dyson (SD) equation with 
the same integral kernel. Our treatment here can include the effects of vector and 
axial-vector mesons. The QCD S parameter is given by the slope of the spin-1 
part of the ''V — A" two-point function, Uwil'^) — ^aa{(i'^), at = 0. We obtain 
the value S = 0.43 ~ 0.48. This is about 30% larger than the experimental value, 
S = 0.31 ~ 0.38 IQ, which comes from the form factors of radiative pion decay 
vr — s> eh"~f and electromagnetic charge radius of pion. Noting that from first Weinberg 
sum rulep the value of this function at = gives the pion decay constant f-„, 
we also calculate the value and show that it is consistent with the previous result 
obtained by solving the homogeneous BS equation for pion.|]TD[| We extract the p 
meson mass and decay constant by three-pole fitting from the two-point function in 
the space-like region (0 < — < (IGeV)^). 



Improved ladder approximation was first used to study the SD equations 11 



12| , and well succeeded to describe the property of the chiral symmetry breaking. 
The homogeneous BS equation in chiral limit was solved in this approximation, 
which represents fundamental properties of pion.||TO|, |13| Especially BS equation 
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for pion was extensively studied in QCD and its generalized model, and the ratio 
between the pion decay constant /tt and the vacuum expectation value {ipip) was 
calculated in Ref. ||lOl- Moreover, the inhomogeneous BS equation in improved 
ladder approximation led to good predictions for low-lying meson masses, fl^ 

We also show the value of S parameter in other dynamical systems. We eval- 
uate S parameter in SU{?>) gauge theory with N£){= 1, ■ ■ ■ , 6) iso-spin doublets of 
fermions in connection with walking technicolor model|]T5[. Although the evaluation 
using dynamical mass function shows that the value of S/No is decreased as No is 



increased []T6[, we find that the value hardly depends on Nd in the improved ladder 
approximation. 

This paper is organized as follows. In section 2, we briefly review the spectral 
representation of S parameter and f.^. We show how to calculate the two-point 
function liyv ~ ^aa from the inhomogeneous BS amphtude. Section 3 is devoted 
to formulations of the inhomogeneous BS equation. In section 4 we show the basic 
tools for the inhomogeneous BS equation and solve it numerically. Section 5 is the 
main part of this paper, where we show the result of the value of QCD S parameter. 
Further we perform three-pole fitting to the two-point function. We also investigate 
the walking coupling case. 



2 QCD S Parameter and Two-Point Function 

In this section, we define the system which we consider in this paper and summa- 
rize the basic ingredients concerning vector and axial- vector two-point functions for 
calculating the QCD S parameter. 

Supposing that u and d quarks are massless in QCD Lagrangian, we have the 
chiral SU{2)l x SU{2)fi symmetry. As is well known, this symmetry is sponta- 
neously broken down to its subgroup SU{2)y, and massless pions appear. Chiral 
Lagrangian well represents the symmetrical aspects of the interaction among pions 
and external currents which couple to photon, W and Z bosons. There are several 
low-energy parameters in the effective chiral Lagrangian and these parameters are 
to be determined from the dynamics of QCD. 

For calculating QCD S parameter we consider the vector and axial- vector current 
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defined by 



T 



u 

d 



(2.1) 



where (a = 1, 2, 3) is Pauli matrix. Then we define the two-point function IIjj 

{J=V,A) 



d^x e^«^ (0|Tj;(x)J^(0) |0) 



(2.2) 



where a, b are iso-spin indices, is the polarization vector defined by e ■ g = 0, 
e - e = -1. 

The QCD S parameter and the pion decay constant are given by the two-point 
function of "V - A" type: 

Air 

dq"^ 

Uvv{0) - IIaa{0) 



S 
f 



Iivy{q^)-IiAA{q^) 



(2.3) 
(2.4) 



QCD S parameter is related to the Gasser-Leutwyler parameter Lio as 5* = — IGvrL 
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Using spectral representation of vector and axial-vector currents, eqs. (^]3|) and 
^41) are rewritten into 



S = A-K — [pAs) - p^{s)] , 

JO 



J TT 



°° ds 



(2.5) 
(2.6) 



where p^(s) and p^(s) are the spin-1 parts of the vector and axial- vector spectral 
functions, respectively. Equations (|2.5| ) and ( p.6|) are referred to as the Das-Mathur- 
Okubo sum rulejl^ and the first Weinberg sum rule||], respectively. We can easily 
understand the above equations by the following way: The commutator of the con- 
served current J^(a;) is decomposed into 



(o| rjx),Jtio) |o) 



A(x: s) 



roo rl Q 

-S'^'l -pAs)isg,, + d,d,)A{x-s) 
d^q 



(2.7) 



-iqx 



(2vr)= 



e{qo)6{s - q 
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where p^^-* denotes a contribution from massless scalar particle. Noting that massless 
pion couples to the axial-vector current ^4° while no massless particle couples to the 
vector current V^", we find p^^^ = fl and p[^^ = 0. The — A" two-point function 
can be expressed as 



(2. 



I J d'x e^''^(0|T [v^{x)V'A^) - Al{x)A'M] |o 

Jo s s — q'^ — le 

Thus, we obtain eqs. ( |2.5| ) and (|2.6|) using eqs. ( p.3|) and ( p.4|) . 
Let us consider the following three-point vertex function 



6i 



2)f J (27r) 



4 e Xa/3ip;q,e) 



■if, 



T^„,^(r/2)^^;(-r/2)j;(a:) 



(2.9) 



where x {=XviXa) bi-spinor, which we call inhomogeneous BS amplitude. Here 
(/) /') ■'O' (^5 h ■ ■ ■) ^-iid (a, (3, ■ ■ ■) denote flavor, color and spinor indices, re- 
spectively. This inhomogeneous BS amplitude has definite spin, parity and charge 
conjugation, i.e., J^*" = 1 for vector case and 1^"'' for axial-vector case. 

Closing the fermion legs of the three-point function, we find that the two-point 
function is expressed in terms of the inhomogeneous BS amplitude x'- 



njj(g^) 



1 



E 



d^p Nc 



tr 



e ■ G) x{p\ q, 



7m 



i(27r)4 2 

for vector vertex 

for axial- vector vertex 



:2.io) 



where Nc = 3 is the number of SU{3)c color and we average over the polarizations for 
convenience; Iljj{q^) in eq. ( |2.2| ) does not depend on the polarization e^. Although 
the vector or axial- vector two-point function itself is logarithmically divergent, the 
chiral symmetry guarantees the finiteness of the quantity nyy(g^) — IlAA{q'^)- The 
divergences appearing in each two-point function cancel because the structures of 



the divergences are exactly the same. |p!8 



The quantities we should first calculate are vector and axial-vector inhomoge- 
neous BS amplitudes Xv Xa > which are finite by current conservation. Then, we 
perform the momentum integration after taking the difference to obtain the — A" 
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two-point function, i.e., 



Ilvv{q^)-IiAA{ 



E 



tr 



:2.1V. 



i(27r)4 2 

This integration converges as we mentioned above. In this paper, we work with 
the space-hke total momentum (g| = —q^ > 0), then we do not encounter the 
singularities which come from meson poles in time-like region. 



3 Inhomogeneous Bethe-Salpeter Equation 

In this section, we discuss the basic formulations for solving the vector and axial- 
vector inhomogeneous BS equations. SD equation is solved by the same BS kernel. 
We give the component form of the inhomogeneous BS equation. The inhomoge- 
neous BS equations are solved in the space-like region for the total momentum g^, 
g|>0. 

In the non-perturbative treatment of QCD by BS approach, the most important 
quantity is the BS kernel which expresses the QCD interaction by the gluon. In the 
improved ladder approximation with Landau gauge, the BS kernel K is defined by 
[the momentum assignment is chosen as shown in Fig. |I| 

Kip, k) = g\p^ ^)Z(^ - ^'^^'^ ^3.1) 

where C2 = (N^ — l)/{2Nc) is the second Casimir of SU{3)c color fundamental 




Figure 1: The Feynman diagram of the BS kernel K{p,k) in eq. ( 3.1 
which we use in our calculation. The helix denotes the gluon propagator 
intermediating between two quarks. 



representation and we use the tensor product notation |[19 

{A ® B)x = AxB . 



(3.2) 
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In the BS kernel ( |3.1|) we adopt the Higashijima-Miransky approximation ||TT|, [T2| 
for the running coupling, g'^{p,k) = 5f^(max(— p^, — A;^)). This running coupling 
allows us to include the property of asymptotic freedom of QCD. Using this running 
coupling the chiral symmetry is always spontaneously broken. The detailed structure 
of the running coupling is given in section [4.1| . 

For solving the inhomogeneous BS equation, we need the full quark propagator 
Sf{p), which is given by solving the SD equation with the same BS kernel for 
preserving the chiral symmetry. The wave function renormalization factor of quark 
propagator is one in Landau gauge in the Higashijima-Miransky approximation. 
Then the quark mass function is given by 

tEip) = K * Sf{p) , (3.3) 

where 

^-(f-) = ray ■ 

The operator "*" acting on the BS kernel K and a bi-spinor denotes momentum 
integration: 

K*^ip)^ l-§^,K(P^k)^{k). (3.5) 
Now, the inhomogeneous BS equation for x{P', <?; e) is 

(T - K*) X = (e ■ G) , ( = 7m, 7m75 ) (3.6) 

where 

T = T{p; q) ^ S]^\p + |) ® S^\p - |) . (3.7) 
The formal solution is given by 

This solution can be reinterpreted by expanding it into the power series of the BS 
kernel as (see Fig.|^) 

X = T-\e ■ G) + T-^K * T-\e ■ G) + T'^ (k * T"^)^ (e ■ G) + • • • . (3.9) 

We can expand the inhomogeneous BS amplitude x = Xv^ Xa i^^*-* eight invariant 
amplitudes: 



X. (p; e) = E (p; g, e) x\ (p; g) , {J = v,A) (3.10) 

i=l 
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Figure 2: The expansion of vector or axial- vector inhomogeneous BS 



amplitude given in eq. (|3.9| ). 



where Xj = 1;''';8) is scalar quantity, r,-'^'' is the vector or axial- vector base 
defined by 

(3.11) 

and 

rl^^ = rr75 , (3.12) 

where = q^j and [a, 6, c] = a [6, c] + 6[c, a] + c[a, h]. We note that the depen- 
dence on the polarization vector is isolated in the base t\^\p\ q, e). 

To solve the integral equation ( |3.6| ) numerically, we perform the Wick rotation on 
the momentum k integration and analytic continuation on the relative momentum 
p as usual. We introduce the scalar variables u, x, v and y as 



p ■ q = —u, = —u"^ — x"^, 
k ■ q = —V, fc^ = — f ^ — y'^ - 



(3.13) 



Multiplying eq. (|3.6|) by the Dirac conjugate base Fj, taking the trace and summing 
over the polarization, we convert it into the component form^ 

J2 {T,j-K,,^)x' = Ii , (3.14) 
i 

where 

Ii{u,x) = j^ti [Ti{p;q,e){e-G) 



^We evaluate the matrix elements of Tij{u,x) and Kij{u,x;v,y) by an algebraic calculation 
program. 
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Tij{u,x) 



7 H [ ri(p; q, e) T(p; q) Tj{p; q, e) 



Kij{u,x;v,y) = 2n J cos 9 -"^ti ^Ti{p;q,e) K{p,k)Tj{k;q,e) ^ . (3.15) 

Here 6 is the angle between the 3-vector part of p and k; cos6 = p ■ k/\p\\k\. The 
operation -k denotes u and x integrations, i.e., 2tc J\ dcosO -k = *. We should note 
that the Dirac conjugation is taken to be 

X{p; q, e) = loXiP*; q\ e)^7o ■ (3.16) 

The complex conjugate on p and q should be taken to preserve Feynman causality 
of inhomogeneous BS amplitudes when those momenta become complex by analytic 
continuation. |jl9|] Our choice of the bases (|3.11|) and ( p.l2|) leads to the fact that the 



matrix Kij is independent of [The q"^ dependence of X'' comes from Tij only.] 

It should be noticed that the invariant amplitude is even function of {p ■ q). 
Namely it is regarded as an even function of u with an arbitrary constant q%: 

X\u,x) = x\-u,x) . (3.17) 

This is the result of the charge conjugation property 

Cx{-P; q, efC-^ = -x{p; q, e) , 
CU-p;q,efC-' = -Up-q,e) , 

where C = z7o72 is charge conjugation matrix. Similarly from this property, one 
can easily check that T and K are real definite and symmetric matrices: 

Tij{u,x) = Tji{u,x) , 
Kij{u,x;v,y) = Kji{v,y;u,x) . (3.19) 

This is an important property, from which we find that the inhomogeneous BS 
amplitudes are real definite, and thus the two-point function is real definite. 



4 Numerical Calculations 

In this section we give the detailed form of the running coupling. We calculate the 
quark mass function and the inhomogeneous BS amplitudes. In the following, all 
the dimensionful parameters are rescaled by Aqcd, otherwise stated. 
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4.1 Running Coupling 

Running coupling can be well approximated by the result from one-loop (3- 

function in high energy region, so we basically use it in the BS kernel (|3.1| ). However 
the one-loop running coupling blows up at /i = Aqcd and we have no idea about the 
functional form in low energy region. One prescription is to adopt the Higashijima 
approximation ini] in which is constant below some scale and is one-loop 



running coupling above that scale. It is important that the running coupling and its 
derivative are continuous with respect to In/i^, otherwise the derivative of the mass 
function is discontinuous.^ We achieve this continuities by interpolating between 
one-loop running coupling and a fixed value with the second order polynomial of 
In /i^. According to Ref. |TD[, we adopt the following functional form of the running 
coupling: 



An 



t 

tp 2tp{tF — tc) 



where t = In /i^ and = 12n/ {llNc — 2Nf) with Nf being the number of flavors. 
As seen in Fig. the support of the — A" two-point function lies below in the 
threshold of c quark, so we take Nf = 3 and ao = An/9. We take the same parameter 
choice as in Ref. [|I^, i.e., tp = 0.5 and tc = —2.0. It is observed at least that the 



specific choice of the parameter tp does not affect the "physical" quantity {ijjijj).flO 
We will check the dependence on the infrared cutoff tp oi the QCD S parameter in 
the section 

4.2 Mass Function 

Before solving the inhomogeneous BS equation, we have to calculate the mass func- 
tion of quarks. The SD equation ( |3.3| ) reads 



= ^ r dy y^^y) (4.2) 

An Jo max(x,y) y-^T?{y) 

in the Higashijima-Miransky approximation, where x = —p^ and y = —k"^. This 

integral equation is solved by the following way: First, we discretize the equation 

^We easily find this point if we convert the SD equation into differential equation. 
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fine enough. Second, starting from tlie functional form S(a;) = constant 7^ 0, we 
iteratively update the mass function according to eq. (|4.2| ) itself until the functional 
form converges. 

When we solve the inhomogeneous BS equation, we need the mass function at 
the momenta p ± q/2. To obtain the mass function S(— ± g/2)^), we substitute 
the above convergent mass function into RHS of eq. ( [4.2|) after putting x = 
— {p ± q/2y, and carry out the y integration. We note that we can independently 
choose the lattice points for solving the SD and inhomogeneous BS equations. 



4.3 Inhomogeneous BS Amplitude 

Let us consider the inhomogeneous BS equation (|3.14|). We discretize it and solve 



the resulting linear equation numerically. Here, we should note that x^i'^iU) is even 
function of v. Then we restrict the integral region over variable v to be positive, 
> 0, after replacing the kernel as 

Kij{u,x;v,y) Kij{u,x;v,y) + Kij{u,x;-v,y) . (4.3) 



The fundamental variables used to solve the inhomogeneous BS equation ( p.l4| ) 
are U and X defined by 

u = exp U , X = expX . (4.4) 

As for the strong interaction, it is important to take into account interactions around 
Aqcd scale rather than that in the high energy scale. The above choice ( [4.4|) is 
suitable for our calculation. We discretize the variables U and X at Nbs = 22 
points evenly spaced in the intervals 

Ue[Xu,Au] = [-5.5, 2.5] , 

(4.5) 

Xe[Ax,Ax] = [-2.5,2.5] . 

In numerical integration, to avoid integrable logarithmic singularity at {u, x) = 
{v,y), we take four-point average of the BS kernel Kij{u,x;v,y) as 



Kij{u,x;v,y) ^[ Kij{u,x;v+,y+) + Kij{u,x;v+,y^) + 

Kij{u,x]v_,y+) + Kij{u,x]v_,y_)] , (4.6) 



10 



where 



v± = exp{V ± -^Du) 



Du = {Au-Xu)/{Nbs-1) , 



(4.7) 



l/± = exp(y ± -Dx) 



Dx = (Ax - \x)KNbs - 1) . 



Now, we solve the inhomogeneous BS equations for the vector and the axial- 
vector currents separately, so that we obtain the amplitudes Xvi'^^^) ^^"^ Xa^'^t^)- 
We use FORTRAN subroutine package for these numerical calculations. 

5 Results 

In this section, first we calculate the spin-1 part of the two-point function, nyv'(g^) — 
nyi^(g^), then extract the QCD S parameter and the pion decay constant /^r. 



After obtaining the vector and axial- vector inhomogeneous BS amplitudes, Xv ^) 
and Xa(^5 numerically, we calculate the "y— A" two-point function using eq. ( p.ll| ). 

The QCD 5* parameter and the pion decay constant are calculated from the 
formulae ( p.3|) and (|2.4| ). Using the numerical differentiation formula, we extract S 
from four data points of livv — ^aa at q\ = —q^ = 0.0, 0.2, 0.4, 0.6. In this choice 
of the interval of q'^, h = 0.2, the error of numerical differentiation is estimated as 
0{h?) ~ 0.8%.|^ On the other hand, the fluctuation 5 of the value of Yiyv — H^a, 
which comes from the dependence on the lattice size, causes the error 0{6/h) of S. 
As we will show below, choosing large lattice size allows us to make the fluctuation 
S within 1%. Then the error of 5* is expected as 0{6/h) ~ 5%. When we take larger 
interval h, the error of numerical differentiation becomes large. On the other hand, 
the fluctuation of S, 0{6/h), is enhanced by taking the smaller interval h. We have 
checked the validity of the differentiation for several choices of h and for various 
numerical differentiation formulae. 

In what follows we investigate all the dependences on the parameters, i.e., the 
momentum cutoff (|4.5|) , the lattice size Nbs and the infrared cutoff tp of the running 
coupling. 

^ It is natural to estimate the numerical error with the dimensionless quantities scaled by Aqcd- 



5.1 QCD 5 Parameter 
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First, we check the dependence on the infrared and ultraviolet cutoff in eq. 
The typical example of the support of the "V — A" two-point function is shown in 
Fig. I with = 0, = 0.5 and Nbs = 22. The choice (gj) covers the dominant 
support very well. Further, the position of the support does not change if we vary 
the values of g^, tp and Nbs- There is a slight rise in the high energy region, which 
originates in the numerical error of the cancellation of divergences between Hyv and 
Haa- This error, if any, affects the two-point function with 1% at most. 




Figure 3: The support of 'V — A" two-point function. The upper 9/10 
of the figure is cripped. 



Second, we check the dependence on the lattice size Nbs- We show the values 
of the QCD 5* parameter and the pion decay constant /^r for several values of Nbs 
in Table |^. We also show the values of Aqcd which are fixed by imposing f^, = 93 
MeV. Even in Nbs = 14 the fluctuation 6 of is within 1% and the fluctuation of 
S is several percents. So, it is enough to take Nbs = 22. 

Third, we check the dependence on the infrared cutoff tp taking A''^^ = 22. 
We show the values of the QCD S parameter and the pion decay constant with 
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Nbs 


10 


12 


14 


16 


18 


20 


22 


S 


0.500 


0.388 


0.452 


0.481 


0.455 


0.474 


0.464 


(/./Aqcd)' X 100 


3.71 


3.76 


3.95 


4.04 


4.04 


4.06 


4.05 


Aqcd [MeV] 


483 


479 


468 


462 


463 


461 


462 



Table 1: Nbs dependence of the value of the QCD S parameter and 
/tt/Aqcd- We also show the values of Aqcd calculated by imposing 
/tt = 93 MeV. We take the parameter choice tp = 0.5. 



several values of tp in Table ^ The variations of S are within 10%. We also calculate 
the scale Aqcd by imposing /^r = 93 MeV, and the results are shown in Table ||. 



The values of Aqcd are consistent with the results in Ref. |T0| which are calculated 
from the homogeneous BS amplitude of pion with the same BS kernel as ours. For 
example, their typical value is Aqcd = 484 MeV with tp = 0.5. [Our tp corresponds 
to tip in Ref. |jTO| as tp = 1 + tjp.] 



tp 


0.3 


0.5 


0.7 


0.9 


1.1 


s 


0.432 


0.464 


0.478 


0.481 


0.470 


(A/Aqcd)' X 100 


3.43 


4.05 


4.10 


3.75 


3.12 


Aqcd [MeV] 


502 


462 


459 


481 


526 



Table 2: tp dependence of the value of the QCD S parameter and 
/tt/Aqcd- We also show the values of Aqcd calculated by imposing 
= 93 MeV. We fix the lattice size as Nbs = 22. 

We show the value of the QCD S parameter which are the main results of this 
paper: 

S = 0.43 ~ 0.48 . (5.1) 

Let us consider what effects are included by our approach. The lowest contribution 
to QCD S parameter is given by one-loop quark diagram shown in Fig. ^(a). There 
are two classes of the higher order corrections to this diagram, i.e., corrections to 
the quark propagator and binding forces to form qq bound state. [Effects of gluon 
condensation [pO|] are not considered here.] These are schematically expressed by the 
diagrams in Figs, ^(b) and H(c). The improved ladder approximation includes these 
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two corrections simultaneously. 




(a) (b) (c) 

Figure 4: The schematic view of the Feynman diagrams contribut- 
ing to the "V — A" two-point function, n(g|;) = IIvv{(1e) ~ ^aa[(^e)- 
(a) is quark one-loop diagram, (b) expresses inclusion of momentum 
dependent mass function, (c) represents binding forces to form quark- 
antiquark bound state. 



QCD S parameters 



Our value 


GL 


VMD 


LCQ 


NCQM 


FQL 


0.43 ~ 0.48 


0.342 ± 0.034 


0.37 


0.21 


0.30 ~ 0.40 


0.16 



Table 3: The predictions of QCD S parameter in various models. GL : 
experimental value by Gasser-Leutwyler, VMD : vector and axial-vector 
dominance model, LGQ : local chiral quark model, NCQM : nonlocal 
constituent-quark model. We also show the value calculated by free 
quark one-loop diagram (FQL), S = Nc/{6Tr). 



Now, let us compare our result with the experimental value (GL)[|I|] and other 

predictions from the vector meson dominance model (VMD) the local chiral quark 

model (LCQ)IP], and the nonlocal constituent- quark model (NCQM)^, |^. We show 

their results in Table The FQL or LCQ model gives a half of the experimental 

^ They calculate I/io(/x) at the renormalization scale fi = rriri (LCQ, NCQM) or fi — nip (VMD). 
Following Ref. 0, |8|, we convert their values of iio(A*) to S by 

^=-16.[w.) + ^(ln^ + l) . 



14 



value (GL) of S. We should include higher order corrections as in Figs, ^(b) and 
§(c). The VMD model well reproduces the experimental value of S. This implies 
that it is important to include the contribution from the bound state, in other words, 
we should take into account the binding force as in Fig. §(c). On the other hand, 
in NCQM model the value of QCD S parameter is improved by the inclusion of the 
momentum dependent mass function as in Fig. Kb). The high energy behavior of 
this diagram is consistent with the result by operator product expansion (OPE).pTl] 
For these reasons we include these two corrections by the improved ladder approxi- 
mation. Our value of QCD S parameter is 30% larger than the experimental value. 
For one thing further corrections beyond the improved ladder approximation may 
be needed (e.g., the decay widths of vector or axial- vector mesons are not taken into 
account in our approximation); for another the difference between our value and the 
experimental one may be caused by the slight breaking of the chiral Ward-Takahashi 
(WT) identity. Reference [p2[ suggests that if we use the ladder approximation com- 
pletely consistent with WT identity we could make up the difference .[] 



5.2 Pole Fitting 

We extract the mass nip and the decay constant fp of p meson by three-pole fitting: 



n(g|) = - 4^ + 4^ , (5.2) 

from the two-point function n(g|;) = liyvi.'fE) ~ ^AA{<fE) in space-like region (0 < 
< (IGeV)^) with tp = 0.5. The positive sign contribution in eq. (|5.2|) comes from 
the two-point function of the vector current and the negative sign contribution from 
that of the axial- vector current. The last two terms represent the contributions from 
other mesons heavier than p meson. The decay constant of p meson in eq. ( ^.2|) is 
defined by 

(0|v;(0)|/(g,e)) = 5%/,m,. (5.3) 

We show our two-point function Il{q%) and the best fitting curve in Fig.|^. The 
best fitted values are shown in Table ^. These values should be compared with those 

^ There is 20% difference in the value of /,r between the "consistent" ladder approximation and 
Pagels-Stokar formula|2^, while Pagels-Stokar formula and the improved ladder approximation 
give almost the same value of /ttQ- [We note that the different BS kernels are used in two 
references. For details see Refs. [0, 
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71. 





(1033 MeV)^ 



Figure 5: Three pole fitting of our two-point function, n(g|;) = 
IIvv{<1e) ~ ^aa{<1e)- The sohd line is the best fitting curve, the dots 
denote the value of our n(g|;). 





Our value 


Experiment [^] 


fp [McV] 


133 


144 ± 8 


rup [MeV] 


643 


768 



Table 4: The best fitted values of the mass and decay constant of p 
meson. We use our value of Aqcd (= 462 [MeV]). 
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in Ref . JTli . Because we need no further regularization in calculating the "V" — A" 



two-point function, we do not have to introduce cutoff parameter as in Ref. ||14 . 

We find that the sum of the pole residues fpTUp — fR^m\_^ + fR^m\^ vanishes, 
which implies that our T^{q\) behaves as I/qe- The result from the improved ladder 
approximation reproduces the high energy behavior of the — A" two-point func- 
tion required by that from OPE in QCD. This means that the spectral functions of 
our n(g|;) in eq. ( |2.8D satisfies the second Weinberg sum rule: 

' ds[pv{s)- pa{s)] = Q . (5.4) 





The masses and decay constants of the heavier mesons are unstable for fitting; 
we obtain several best fitted curves with different values for the masses and decay 
constants of the heavier mesons, although they satisfy the first and second Weinberg 
sum rules. On the other hand, the lowest meson mass and decay constant are very 
stable. 

We compare our two-point function with that from the p and ai meson dominance 
model. In this model the first and second Weinberg sum rules read 

fl - fl = fl , - fK = . (5.5) 

It is convenient to adopt the following parameterization: 

fp = /tt cosh e , fai = sinh 6 . (5.6) 
Then the resultant "V^ — A" two-point function F^q]^) is expressed by 

. 2. 2 m^coth^^ ^ ^ 

(g| + mj) [q% + coth 0) 

We show the function F{ql^) for two different choices of mp and coth 6' in Fig. |^ with 
our n(g|;). Taking (coth 6', m^) = (1,643 [MeV]) makes Flq"^) well agree with our 
n(g|;). However, when we use the experimental value {fp,mp) = (144, 768) [MeV], 
F[q^) does not match with our Ii{q\). 

5.3 Walking Coupling Case 

We apply the method for inhomogeneous BS equation to the other dynamical system 
than QCD, i.e., SU{?>) gauge theory with A^^, doublets of massless fermions. 
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Figure 6: The comparison with the p and ai meson dominance modeL 
The sohd hne denotes F{q%) with (coth^,mp) = (1,643 [MeV]), the 
broken hne, with {fp.mp) = (144,768) [MeV]. The dots denote 

the values of our n(g|;). 
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The difference between the real QCD and the system which we investigate in this 
section appears in the coefficient of the /^-function. Namely, the coefficient factor 
tto of the running coupling in eq. ( [4.1|) is taken to be 

127r 



An 
To 



llNr - AN, 



(5i 



D 



where Po is the coefficient of one-loop /?-f unction. For the small value of Po the 
coupling runs very slowly, and the system well simulates the walking technicolor 
model||15|. 

Using the same procedure as the previous one for real QCD, we evaluate S 
parameter and the pion decay constant with N^, = 1, ■ ■ ■ , 6. The results are shown 
in Table |^. We find that the dominant support (cf. Fig. |^) lies on the energy region 





1 


2 


3 


4 


5 


6 


S/Nd 


0.469 


0.457 


0.442 


0.431 


0.430 


0.421 


(/JAqcd)VA^d X 100 


3.87 


4.27 


4.87 


6.05 


9.32 


24.7 


Aqcd [MeV] 


473 


318 


243 


189 


136 


76.4 



Table 5: The values of S/Nd and /tt/Aqcd for the number of doublets 
Nd = 1, ■ ■ ■ , 6 with tp = 0.5 fixed. We also show the values of Aqcd by 
imposing = 93 MeV. 

higher for large N^ than for small Nd. The values of Aqcd become small for large 
number of doublets. However, we cannot see the clear dependence on Nd of the 
values of S/Nd- The value slightly decreases as Nd increased, only by 11%. Our 
result shows that S parameter almost scales linearly with Nd, which supports the 
naive estimation by Refs. 0, 

Finally we check the dependence on the infrared cutoff tp with Nd = 5, and 
show the results in Table |^. The dependence on the infrared cutoff tp of Aqcd 



should be compared with the result in the case of single-sextet quark in Ref. [|T0 
They calculate /tt/Aqcd using the homogeneous BS equation of pion, and find that 
the slowly running coupling gives stable results less dependent on tp. Our results 
conform to theirs. 
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tp 


0.3 


0.5 


0.7 


0.9 


1.1 


S/Nd 


0.450 


0.430 


0.432 


0.436 


0.446 


(A/Aqcd)V^d X 100 


8.82 


9.32 


9.78 


10.2 


10.4 


Aqcd [MeV] 


140 


136 


133 


130 


129 



Table 6: The dependence on the infrared cutoff tp for Nd = 5. 
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